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Abstract 

We present detailed results of a diagrammatic calculation of the leading two-loop 
QCD corrections to the masses of the neutral CP-even Higgs bosons in the Minimal 
Supersymmetric Standard Model (MSSM). The two-loop corrections are incorporated 
into the full diagrammatic one-loop result and supplemented with refinement terms 
that take into account leading electroweak two-loop and higher-order QCD contribu- 
tions. The dependence of the results for the Higgs-boson masses on the various MSSM 
parameters is analyzed in detail, with a particular focus on the part of the parameter 
space accessible at LEP2 and the upgraded Tevatron. For the mass of the lightest 
Higgs boson, m^, a parameter scan has been performed, yielding an upper limit on 
which depends only on tan/3. The results for the Higgs-boson masses are compared 
with results obtained by renormalization group methods. Good agreement is found in 
the case of vanishing mixing in the scalar quark sector, while sizable deviations occur 
if squark mixing is taken into account. 
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1 Introduction 



The search for the lightest Higgs boson is a crucial test of Supersymmetry (SUSY) which 
can be performed with the present and the next generation of accelerators. The prediction 
of a relatively light Higgs boson is common to all supersymmetric models whose couplings 
remain in the perturbative regime up to a very high energy scale [[[]. A precise prediction for 
the mass of the lightest Higgs boson in terms of the relevant SUSY parameters is necessary 
in order to determine the discovery and exclusion potential of LEP2 and the upgraded 
Tevatron, and for physics at the LHC and future linear colliders, where eventually a high- 
precision measurement of the mass of this particle might be possible. A precise knowledge 
of the mass of the heavier CP-even Higgs boson, m#, is important for resolving the mass 
splitting between the CP-even and -odd Higgs-boson masses. 

In the Minimal Supersymmetric Standard Model (MSSM) M at the tree level the mass 
m/j of the lightest Higgs boson is restricted to be smaller than the Z-boson mass. How- 
ever, this bound is strongly affected by the inclusion of radiative corrections: the dominant 
one-loop corrections arise from the top and scalar-top loops which yield terms of the form 
G Firif Inim^m^ / m%) |§. These results have been improved by performing a complete one- 
loop calculation in the on-shell scheme, which takes into account the contributions of all sec- 
tors of the MSSM ||, j|, ||. Beyond one-loop order, renormalization group (RG) methods have 
been applied in order to include leading logarithmic higher-order contributions |7|, |8|, ||, [UJ . 
In the effective potential approach diagrammatic results for the dominant two-loop contri- 
butions have been obtained in the limiting case of vanishing t-mixing and infinitely large 
Ma and tan/5 |IT |. The calculation of the leading QCD corrections in this approach has 



recently been generalized to the case of arbitrary tan/? and non- vanishing t- mixing |T^|. 

Up to now phenomenological analyses have been based either on the RG results [0, ||, 
[| [lO], or on the complete one- loop on-shell results @, ^ §]. These results differ by large 
leading logarithmic higher-order contributions, which are not included in the one-loop on- 
shell results, but also by non-leading one-loop contributions, which are neglected in the 
RG approach. The numerical difference in the Higgs-mass predictions between the two 
approaches reaches up to 20 GeV. 

Recently a Feynman-diagrammatic calculation of the leading two-loop corrections of 
0(aa s ) to the masses of the neutral CP-even Higgs bosons has been performed [13|, [14 . 



Compared to the leading one- loop result the two-loop contribution was found to give rise 
to a considerable reduction of the value. The leading two-loop corrections have been 
combined with the full diagrammatic one-loop on-shell result || and further refinements 
have been included concerning the leading two-loop Yukawa corrections of 0(G 2 F mf) || pl| 
and leading QCD corrections beyond two-loop order. 

In this paper we present in detail the steps of this calculation. The results for the masses 
of the neutral CP-even Higgs bosons are analyzed in terms of the relevant parameters of the 
MSSM. A parameter scan for the lightest Higgs-boson mass is performed yielding an upper 
bound for uih within the MSSM (apart from certain threshold regions which correspond to 
very specific configurations of the MSSM parameters) given exclusively in terms of tan/3 . 
This upper bound is discussed in view of the discovery potential of LEP2 and the upgraded 



1 



Tevatron. The results for are compared with the corresponding results obtained by RG 
methods. The comparison is performed both in terms of the (unobservable) parameters of 
the scalar top mass matrix and in terms of the physical stop masses and the stop mixing 
angle. 

The paper is organized as follows: Section 2 contains our notations and a description 
fo the renormalization procedure as required for the corrections in the MSSM Higgs sector 
in 0(aa s ). The main features of the calculation are discussed in section 3. In section 4 
we present a detailed numerical analysis of the results for the neutral CP-even Higgs-boson 
masses as functions of the different SUSY parameters. We perform a scan for over the 
parameters nig, Ma, M, li and the t-mixing parameter and determine the maximal possible 
values of as a function of tan (3 . Finally numerical comparisons are shown with results 
obtained by renormalization group (RG) methods. In section 5 we give our conclusions. 



2 Renormalization 

2.1 The Higgs sector of the MSSM 

The Higgs sector of the MSSM consists of two Higgs doublets Hi, H 2 with opposite hyper- 
charges Yi — — 1 and Y 2 = +1 and non- vanishing vacuum expectation values V\ and v 2 . The 
Higgs doublets can be decomposed according to 



Hi 
Ho 



H°i ^ = (vi + ^{cj>l + i X ^ 

h)\ = ( 

Hi ) ~ \v 2 + ±{^ + ixl) 



V2^ 



The vacuum expectation values define the angle j3 via 



Vo 

tan/3 = — ; < (3 < it/2 . (2) 

Vl 



The Higgs potential, including all soft SUSY breaking terms reads |L6| (ei 2 = — 1): 

V = m\\Hi\ 2 + m 2 2 \H 2 \ 2 -m 2 l2 (e ab HlH b 2 + \i.c.) 

+ \{9i 2 + 92 2 ) [m 2 - \H 2 \ 2 } 2 + l -g 2 2 \R\E 2 \ 2 , (3) 

where m 2 = \n\ 2 +fh 2 {i = 1,2); rhi, rh 2 , m i2 are the soft SUSY breaking terms, and \x denotes 
the mixing between Hi and H 2 . The coupling constants of the Higgs self- interaction are, 
contrary to the SM, determined through the gauge coupling constants gi and g 2 . Besides 
gi, g 2 two independent parameters are required to fix the potential (0) at the tree level. 
Conventionally they are chosen as tan/3 and M 2 A = — m^ 2 (tan/3 + cot f3 ), where Ma is the 
mass of the CP-odd A boson. 
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The diagonalization of the bilinear part of the Higgs potential, i.e. the Higgs 
trices, is performed via the orthogonal transformations 



mass ma- 



h° 
' G° 

, A ° 
G ± 



cos a 
— sin a 

cos ft 



sma 
cos a 
sin ft 



— sin ft cos ft 
cos ft sin (3 

— sin (3 cos ft 



0? 

xS 

4>t 



with ft from eq. (0). The mixing angle a is determined through 

M\ + M 2 Z _ 



tan 2a = tan 2ft 



71 



< a < 



One gets the following Higgs spectrum: 



2 neutral bosons, CV = +1 
1 neutral boson, CV = — 1 
2 charged bosons 
3 unphysical Goldstone bosons 



h° } H° 
A 

H + ,H- 

r 1 ' 



The masses of the gauge bosons are given in analogy to the SM: 



Mi 



u = 2^2(^1+^2); 



Mi 



2^1 +fe 2 )K 2 + ^ 2 ); 



(4) 
(5) 
(6) 



(7) 



Mv = 0. 



(9) 



At tree level the mass matrix of the neutral CP-even Higgs bosons is given in the 0i- 
basis in terms of M z , Ma, and tan/3 by 



M 



2, tree 
Higgs 



™>4 
ml 



m- 4 
m 



2 



Mi sin 2 /? +M|cos 2 /3 
-(Mi + Mf) sin/? cos/? 



-(Mi + Ml) sin/? cos/? 
Mi cos 2 ft + Ml sin 2 /3 



(10) 



which by diagonalization according to eq. (f|) yields the tree-level Higgs-boson masses 



1V1 Higgs 



if, tree 







Til 



m 



/i.trec 



In order to slightly simplify the two-loop calculation, we have chosen to perform it in 
the 0i-02-basis. In this way the angle a does not appear in the calculation of the two-loop 
self-energies, but enters at the end when the rotation into the physical basis is performed. 
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In order to deal with the arising divergencies and to establish the meaning of the phys- 
ical parameters beyond the tree level, one has to renormalize the Higgs and the scalar top 
sector of the MSSM. For the corrections of 0(aa s ) to the Higgs-boson masses, in the focus 
of this discussion here, renormalization up to the two-loop level is needed. In the following 
we specify the renormalization for the relevant quantities in this calculation (explicitly listed 
are only those terms that actually contribute at 0(aa s )). The renormalization of the com- 
plete one-loop contributions to the neutral CP-even Higgs-boson masses has been performed 
according to Ref. 0. 

We use the following notation: E^ and E^ 2 ) denote the one- and two-loop part of an 
unrenormalized self-energy, EW and E^ 2 ) denote the one- and two-loop part of a renormalized 
self-energy, and E'(A; 2 ) = t^E^ 2 ) . t\ and t 2 denote the unrenormalized tadpoles; and 
£( 2 ) represent the one- and two-loop part of an unrenormalized tadpole, and ti,t 2 denote the 
renormalized tadpoles. 

The renormalization of the masses and fields is performed as follows: 



Mf M 2 Z + 5M 2 Z {1) + 5M 2 Z {2 \ (12) 

Ml M 2 A + 5M 2 A (1) +5M 2 A {2 \ (13) 

ipi -> <PiZ% 2 , </?i = 0?,x?,0r> ( 14 ) 

(f 2 -> <f2Z]l 2 , ip 2 = (f>2,X2, <f>t> ( 15 ) 

Z Hi = 1 + bZl] + 5Z%] (16) 

tan/3 -> tan/3(l + 5tan/3 (1) +5tan/3 (2) ), (17) 

U -> U + 5t\ 1} + 5t {2) (i = 1,2). (18) 

This yields for the renormalized two- loop self-energies of 0i and <p 2 '- 

tf,{k 2 ) = ^(k 2 ) + k 2 5Z$-5v£\ (19) 

Sg(* 2 ) = X%(k 2 ) + k 2 5Z%l-8v£\ (20) 

S^ 2 (^ 2 ) = 41(^ 2 )-< 2 L (21) 



where it is understood that the unrenormalized self-energies at two-loop order also contain 

(2) (2) 

the contributions arising from the subloop renormalization. The expressions o , oV<L and 
SVli are the two-loop counterterm contributions from the Higgs potential: 

SV} 2) = 5M 2 Z {2) cos 2 (3 + 5M 2 A {2) sin 2 p - 5t {2) cos (3 (1 + sin 2 (3 ) 

v 2M w sw 

+^2 ] 777~f~ cos 2 p sin p + 5Z { hI (m 2 z cos 2 p + Ml sin 2 p ) 

+5 tan p (2) cos 2 P sin 2 p{M 2 A -M 2 z ), (22) 
5Vj 2) = 5M 2 Z (2) sin 2 p + 5M 2 A {2) cos 2 p - 5t { 2 2) — ^ sin p (1 + cos 2 p ) 
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+5t {2) — ^ sin 2 p cos p + 6Z& (M 2 z sin 2 p + M 2 cos 2 p ) 

-<5 tan /? (2) cos 2 p sin 2 /3 ( Af| - M|) , (23) 
tfV^ = -5M 2 Z {2) sin /? cos /? - 5M 2 A {2) sin /? cos p - 5tf ) - sin 3 p 



2M w s w 

-5tf C cos 3 /3 - *zg Sin/j C ° S(3 (M\ - M 2 Z 
2 2M w s w H Hl 2 V 



2 

-5 tan /3 (2) 1 sin 2/5 cos 2P [M 2 A + M|), (24) 
with the electroweak mixing angle s^y = 1 — c 2 ^, c 2 ^ = M^/M|. 

The counterterms are fixed by imposing on-shell renormalization conditions for the renor- 
malized self-energies. For the A boson this reads: 

Ret A (Ml)=0. (25) 

The tadpole conditions are: 

*i = 0, t 2 = . (26) 

The conditions for the tadpoles have the consequence that the t>j remain the minima of the 
Higgs potential also at the two-loop level. 

The resulting expressions for the renormalization constants contributing to the leading 
two-loop corrections to the neutral CP-even Higgs-boson masses, expressed in terms of un- 



renormalized self-energies and tadpoles, are given in Sec. 3.1 



2.2 The scalar quark sector of the MSSM 

Renormalization in the squark sector is needed in the present calculation at the one-loop 
level, i.e. at 0(a s ). As above, we work in the on-shell scheme. In the following the formulas 
are written for one flavor. 

The squark mass term of the MSSM Lagrangian is given by 

C mi = -\(flP R )z ( { L ] , (27) 



where 



MJ + Mf cos 2P (l( - Q f s 2 w ) + mj m f (A f - //{cot P ;tan/3 }) 



Z= | Q * 1 v ° ^ 1 J " 1 1 L ' ' | • (28) 

m f (A f - //{cot /3; tan /3}) Mj, + Mf cos 2/3 Q f s 2 w + mj 

and {cot/3 ;tan/3 } corresponds to {u; <i}-type squarks. The soft SUSY breaking term Mq, 
is given by: 

„ J Mjj for right handed u-type squarks , , 

Ql 1 Mf) for right handed <i-type squarks 



In order to diagonalize the mass matrix and to determine the physical mass eigenstates 
the following rotation has to be performed: 



fi 

h 



cos Qj sin Of 
— sin Oj cos Oj 



fn 



The mixing angle 8? is given for tan (3 > 1 by: 



(30) 



cos 9? 



sin Of 



\ m 2 (Af — //{cot (3 ; tan/3 }) 2 + (m 2 ^ — m 2 ^) 2 



(31) 



=F sgn 



Af — //{cot (3 ; tan /? } 



x 



m 2 (Af — //{cot /3 ; tan/3 }) ; 



\ m) (A/-//{cot/?;tan/?}) 2 + 



m 2 ? — m 2 ? ) 2 

j J? ,/ 1 



(32) 



The negative sign in ([52]) corresponds to u-type squarks, the positive sign to <i-type ones. 
M? + M| cos 2/3 Q fSyy + m 2 denotes the lower right entry in the squark mass 



matrix (|28|) . The masses are given by the eigenvalues of the mass matrix: 



m , 



Ml + M% 



+ -Mf cos 2/3// + m 2 f 



± 








± 









c f 



sgn 



M? - MJ, + Mf cos 2/3 (J| - 2Q / s 2 v )] + 4m 2 (A u - // cot /3 

Mj - Mj, + M| cos 2/3 (// - 2Q / s 2 y )] 2 + 4m 2 (A d - //tan/3 
MJ - Mj, + Mf cos 2/3 (// - 2Q f s 2 w 



.(33) 



(34) 



for w-type and <i-type squarks, respectively. For most of our discussions (see Sec. £|) we make 
the choice 



M 



Q 



Mr,, =: m. 



>j ■ 



(35) 

Since the non-diagonal entry of the mass matrix eq. ( p8|) is proportional to the fermion mass, 
mixing becomes particularly important for / = t , in the case of tan f3 ^> 1 also for / = b. 

For an on-shell renormalization it is convenient to express the squark mass matrix in 
terms of the physical masses mj^m^ and the mixing angle 0j: 



cos 2 6fm% + sin 2 6fm% sin 6? cos Offm 2 ; — m 2 
J fi j h j j v ji 

sin Of cos Of {m 2 j — m 2 J sin 2 ^m 2 +cos 2 ^m 2 



(36) 



7"7i 



"7 "7a 



can be written as follows: 



Ai 



sin cos 0f(m 2 7 — m\ ) 
J j v ji HJ_ 

m f 



+ //{cot (3 ; tan /3 }. 



(37) 



The renormalization of the fields, the masses, and the mixing angle is then performed 



via 



h 

2 



m 



0, 



7 



-SZf ) 
mf + 5m 2 ? 

Ji Ji 



In the mass eigenstate basis, the field renormalization reads: 



fi 
h 



1 + \5Z h 



1 + ¥ z h 



h 
h 



with 



5Z 



5Z 



h 



h 
5Z 



cos 2 9j sin 2 Qj 
sin 2 9j cos 2 0j 



5Z 
SZ 



h 

Jr 



fl2 



= sin Oj cos 0j(8Zj r — 6Zf L ) = 5Zf 2 _ 



(38) 

(39) 

(40) 
(41) 



(42) 

(43) 
(44) 



sin 0j cos Of 
cos 2 Of — sin 2 Of 



-(5Zf 2 -5Z h ). 



The renormalized diagonal and non-diagonal self-energies in this basis have the following 
structure: 



± h (k 2 ) = Y, h {k 2 )-5ml + {k 2 -m 2 h )5Z h 
t h (k 2 ) = H h {k 2 )-5m 2 h + {k 2 -m) 2 )5Z h 

Hh^) = %/ 2 (* 2 ) - K - ™%)88f + (* 2 - \{m 2 h + m\))5Z h 
We impose the following on-shell renormalization conditions: 
Re XL 



(45) 
(46) 

(47) 



Re t' f 
ReE /; 
ReE'~ 



ReE f ~ f 









( m 2 \ = _ Re E / / m 2 \ + Re S / / m 2 \ 

v Jl ' Jl y Jl ' J2 K S2> 

K) = o, 

which determines the renormalization constants to be 



5m ; = ReEf (m,-) 
5m 2 ? = Re E ? (m| ) 

72 /l V 72 ' 



(49) 
(50) 
(51) 
(52) 

(53) 
(54) 



sz } 
sz 



-E'r (m% 



h 
56? 



8Z f 



6Z f 

112 



m ; — m • 
fi h 



(55) 
(56) 

(57) 



The unsymmetric renormalization condition (|5l|) is chosen for convenience since it leads to 
5Zj 2 = SZ^ and accordingly to 5Zi = 0, which simplifies the expression for the counterterm 



of the mixing angle. In eq. fl52|) we have imposed the condition that the non-diagonal self- 



energy vanishes at q 



ml 



Alternatively one could choose q 



m z , instead; the numerical 

H 



difference arising from these different choices is irrelevant for the results of the Higgs-boson 
masses, as we have checked explicitly. 

Taking into account that neither 5fi nor 5 tan/? are of 0(a s ), one obtains from eq. 



5 A, 



sin 9? cos 6?(m 2 7 — mh 
j j h /2 

m f 



1 - 2 sin 2 e~ f 5m 2 - f 

f se f - + — £- 



5m 2 ; 
h 



sin Qj cos 9j 



m 



h 



m 



h 



5m 4 



m\ 



(5? 



For completeness we also list the expression for the quark mass counterterm in the on- 
shell scheme, 

tfm / = m / (E]T(TOj) + Ej(m})) l (59) 

where the scalar functions in the decomposition of the fermion self-energy E/(p) are defined 
according to 

E f (p) =^(p 2 ) +^f(p 2 ) +m / Sj(p 2 ). (60) 



3 Calculation of the neutral CP-even Higgs-boson 

masses 

3.1 Leading two- loop contributions to the Higgs-boson 
self-energies 

The dominant one-loop contributions to the Higgs-boson mass matrix in eq. ( |10"1) are given 
by terms of the form Gpm\ In (m^m^/m 2 ), which arise from t- and t-loops. They can be 
obtained by evaluating the contribution of the t-t-sector to the self-energies at zero 
external momentum from the Yukawa part of the theory (neglecting the gauge couplings). 
Accordingly, the leading contributions to the one-loop corrected Higgs-boson masses are 
derived by diagonalizing the matrix 

M 2,i-ioo P _ m l,~ £^(0) ml 1<h - S^ 2 (0) \ ( . 

where the E^ denote the one-loop Yukawa contributions of the t-t-sector to the renormahzed 
one-loop 0x,2 self-energies. For completeness, we list here the explicit form of these dominant 
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one-loop corrections (in the numerical results given in Sec. ^ we use the complete one-loop 
on- shell result as given in Ref. ||): 



,(i) 



,(o) 



4>) 



H m |) 2 V 

fi(A t — fi cot (3 



SGpmf [i 2 (A t — fi cot (3 
V2n 2 sin 2 (3 (m 

2y/2n 2 sin 2 /? »ir 



2fj,A t (A t — \i cot /3 
(m|-m|) 2 



2 

m~ 

12 



2 

ti 

2 

mi 
In— ^ 



2 
*2 



1 - 



2 2 

— mf 

*1 *2 



m *2. 



hi 



V2n 2 sin 2 /? 

| ^ 2 (A -/icot/3 
(m| - m|) 2 



1 



| A f (A t -/fCOt/3) 

/ m 2 — m~ 

*1 *2 

m l + m l ln ^? 1 
m~ — m~ 

tl 52 



hi 



2 
El 



m 



t2 



m 



(62) 



By comparison with the full one-loop result 0, [5], [(J it has been shown that these con- 
tributions indeed contain the bulk of the one-loop corrections. They typically approximate 
the full one-loop result within 5 GeV. 



In order to derive the leading two-loop contributions to the masses of the neutral CV- 

13|, inj. Accordingly, 



even Higgs bosons we have evaluated the QCD corrections to eq. (pip 
we have calculated the 0(aa s ) contribution of the t-t-sector to the 0^2 self-energies at zero 
momentum transfer, neglecting the gauge couplings. Because of the large value of the strong 
coupling constant these are expected to be the most sizable two-loop corrections (see also 
Ref. 



The leading two-loop contributions to the 01,2 self-energies are given, according to 



Eqs. 



ID, by 



4 2 >) 




(63) 


sg(o) 


= sg(o)-< 2) , 


(64) 






(65) 



and for the leading contributions the potential counterterms eqs. (22)-(^4|) simplify to 



sv 



(2) 



SV 



(2) 



sv 



(2) 



+5M 2 A (2) sin 2 (3 -St\ 



(2) 



2M w sw 



2M w s w 



cos 2 (3 sin (3 , 



+5M 



2{2) cos 2 /3 -M 2) 



2M w sw 



cos (3 (1 + sin 2 (3 



sin (3 (1 + cos 2 (3 



+St\ 



(2) 



2M w sw 

r 2 (2) 



sin 2 (3 cos (3 



SMX { > sin 13 cos/3 - St\ 



(2) 



2M w s w 



sin 3 /3 



(66) 



(67) 



(6* 



-8tf— ^ cos 3 /3. 

2M w sw 

From the on-shell renormalization conditions eqs. (f25|)-(|26|) we obtain for the countert- 
erms in eqs. (|BB|)-(|5SD 

5M 2 A {2) = £^(0) (69) 

and 

M 2) = -4 2) , M 2) = -4 2) - (to) 



3.2 Evaluation of the relevant Feynman diagrams 

The calculations have been performed using Dimensional Reduction (DRED) W7\, which is 
necessary in order to preserve the relevant SUSY relations. Naive application (without an 
appropriate shift in the couplings) of Dimensional Regularization (DREG) [|18| , on the other 
hand, does not lead to a finite result. The same observation has also been made in Ref. JTl] . 



The Feynman diagrams contributing to the 0i,02 and A self-energies are depicted in 
Fig. |l].[] The Feynman diagrams for the tadpole diagrams are shown in Fig. |2]. 

There are three classes of diagrams: pure scalar diagrams (Fig. [TJa-c, Fig [2]a), diagrams 
with gluon exchange (Fig. [l]d-h, Fig |2|b-c), and diagrams with gluino exchange (Fig. [p-l, 
Fig |]d-e). These diagrams have to be supplemented by the corresponding one-loop diagrams 
with counterterm insertions, which are depicted in Fig. |3| and in Fig. f|. The counterterm in- 



sertions are generated by the renormalization in the top and scalar top sector (see Sect. 2.2). 
They are calculated from the Feynman diagrams in Fig. |5]. 

The gluon-exchange contribution of 0(a s ) to the quark mass counterterm reads in DRED: 

69mq = ^l mq fj_ + 7E + l n (J!!!L) (S), (71) 



where 25 = 4 — n with n the space-time dimension, 7e is Euler's constant, and /i is the 
't Hooft scale. The explicit form of the other counterterms of the quark and scalar quark 
sector can be found in Ref. |TP| . 



Some of the diagrams shown in Figs. [I], @ vanish when they are combined with the 
corresponding counterterm contributions of Figs. |3|, [|. From the pure scalar diagrams only 
Fig. |l]a yields a non-vanishing contribution. The diagrams Fig. [l]b-c are canceled exactly 
with their corresponding counterterm diagrams. Here the mass renormalization for the 
diagonal terms (with two identical squarks) and the mixing-angle renormalization for the 
non-diagonal terms (with two different squarks) are needed. The same applies for the tadpole 
diagram Fig. together with the counterterm diagram Fig. |]b. The diagrams Fig. [Ijf are 
exactly canceled with the corresponding diagram with counterterm insertion Fig. ^|b. The 
same applies for the tadpole diagrams Fig. [2]b together with the counterterm diagram Fig. |]b. 

1 The diagrams with a closed gluon line give zero contribution in DREG and DRED, they are omitted 
here. 



10 



We now briefly describe the evaluation of the two-loop diagrams. As explained above, 
the calculation involves irreducible two-loop diagrams at zero momentum-transfer and coun- 
terterm diagrams. In deriving our results we have made strong use of computer algebra 
tools: the diagrams were generated with the Mathematica package FeynArts [pOfl . F° r this 
purpose we have implemented a model file which contains the relevant part of the MSSM 
Lagrangian, i.e. all SUSY propagators (ti, t^, bi, &2, g) needed for the QCD-corrections and 
the appropriate vertices (Higgs boson-squark vertices, squark-gluon and squark-gluino ver- 
tices). The program inserts propagators and vertices into the graphs in all possible ways 
and creates the amplitudes including all symmetry factors. The evaluation of the two-loop 
diagrams and counterterms was performed with the Mathematica package TwoCalc |2T|. By 
means of two-loop tensor integral decompositions it reduces the amplitudes to a minimal set 
of standard scalar integrals, consisting in this case of products of the basic one-loop integrals 
Aq, B (the B functions originate from the counterterm contributions only) and the 



two- loop function T134, which is the genuine two-loop scalar integral at zero momentum- 
transfer (vacuum integral). This integral is known for arbitrary internal masses and admits 
a compact representation for 5 — > in terms of logarithms and dilogarithms (see for instance 
p3||). It should be noted that from the expansion of the one-loop two-point function B , 

1 



Ref. 



B (q ,m a ,m b ) 



Bl n (q\ m a , m b ) + 5 B 6 (q\ m a , m b ) 



5 



(72) 



only the term B$ n contributes, while Bq drops out in our final result. From the output 
generated with TwoCalc a FORTRAN code was created which allows a fast calculation for 
a given set of parameters. This code has been implemented into the FORTRAN program 
FeynHiggs fl2"4jl , see below. 

Our results for the two-loop <pi,2 self-energies are given in terms of the SUSY parameters 
tan/3 , M4, ji, m^, m f ~ 2 , 6* f -, and m g . In the general case the results are by far too lengthy to 
be given here explicitly. In the special case of vanishing mixing in the t-sector, /i = 0, and 



m 



m 



t-2 



rrif, a relatively compact expression can be derived which is given in Ref. 



We have performed an expansion of this result for large values of trig. It yields for the leading 
terms 



,(2) 

J <l>2 



a a. 



3 mr 



71 71 2 Sw 



sin 



-In 2 



+ ln(m?) +C(m? 



(73) 



This shows that the gluino does not decouple from the two-loop result, contrary to the case 
of the two-loop QCD contributions to the p-parameter in the MSSM |25[| . 
In Ref. JTT] a result for the limiting case 

mi = = mj 2 = 3> m t , tan/? — > 00 (74) 

has been given. In this limit we obtain 



4 2 >) = °> Egfc(O) 



(2) 



4?(o) 



-a 



a a, 



9 m\ 



7T 7T AS 2 W M^ 



111 



mi 



In [ + 2 



m 



(75) 



which agrees with the corresponding result given in Ref. [IT|. 
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Figure 1: Feynman diagrams for the contribution of squark and quark loops to the Hig 
boson self-energies at the two-loop level (H = 0i, 2 , A). 
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Figure 2: Feynman diagrams for the contributions of squark and quark loops to the Higgs- 
boson tadpoles at the two-loop level (H — <f>i, 2 ). 
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Figure 3: One-loop counterterm contributions to the Higgs boson self-energies at the two- 
loop level (H = <f>i,<f> 2 ,A). 
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Figure 4: One-loop counterterm contributions to the Higgs-boson tadpoles at the two-loop 
level (H — 0i, 2 )- 
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Figure 5: One-loop diagrams for the squark and quark mass counterterms and for the 
squark mixing-angle counterterm. 
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3.3 Determination of the Higgs-boson masses 



In the Feynman-diagrammatic approach the Higgs-boson masses are derived beyond tree 
level by determining the poles of the h — //-propagator matrix whose inverse is given by 



(A 



HiggsJ 



2 2 

q — m HtTee 



Z hH (q 2 ) 



Z hH (q 2 ) q 2 ~ ml tree + Z h (q 2 



(76) 



where again the 2 denote the renormalized Higgs-boson self-energies, now in the h — H -basis. 
Determining the poles of the matrix AHi ggs in eq. (ff6|) is equivalent to solving the equation 



q 2 - m\ + t hh (q 2 ) q 2 -m 2 H + t HH (q 2 



%iu{q 2 



0. 



(77) 



In our calculation the complete one-loop result for the Higgs-boson self-energies in the 
on-shell scheme || is combined with the leading two-loop contributions, which have been 
outlined in the previous section. The matrix eq. ([76]) therefore contains the renormalized 
Higgs-boson self-energies 



£ s (q 2 ) = £ ( P(q 2 ) + ±?\0), s = h,H,hH, 



(7E 



where the momentum dependence is neglected only in the two-loop contribution. 

Since the two-loop contribution has been calculated in the </>i-02-basis, a rotation into 
the /i- if -basis, according to eq. (01), has to be performed: 



M2) 
-'H 

M2) 

-'h 

,(2) 
J hH 



2 v(2) , • 2 v(2) , o • (2) 

cos a + sin a L^ 2 ' + 2 sm a cos a 

• 2 v(2) . 2 v(2) o • v(2) 

sm a L\ / + cos a LV — 2 sm a cos a LA/, 



— sm a cos a 



(*8? - s 



+ (cos a — sin a )£ 



(2) 



(79) 



We have implemented two further corrections beyond 0(aa s ) into the prediction for 
m h , which are illustrated in Figs. §, [7|, |9] and [!(]. The leading two- loop Yukawa correction 
of 0(G 2 F mf) is taken over from the result obtained by renormalization group methods. It 
reads |H 111 



Am 2 h 



with 



9 



167T 

X 



,G 2 F m\ 



xt + r 



1 2 

mf — mf 

t2 I 



Am 2 



+ 



9 

mf 

52 



log 



2m 2 
'm 2 ^m 2 2 
mf 



3- sin 2 20* 



sin 2^- log 




30) 



(81) 
(82) 
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The second higher-order contribution which has been implemented concerns leading QCD 
corrections beyond two-loop order, taken into account by using the MS top mass0 



1 + j^Usijnt) 



m t = mt{m t ) « - ; 4 — N (83) 



for the two-loop contributions instead of the pole mass, mt = 175 GeV. In the t mass matrix, 
however, we continue to use the pole mass as an input parameter. Only when performing 
the comparison with the RG results we use m t in the t mass matrix for the two- loop result, 
since in the RG results the running masses appear everywhere. This three-loop effect gives 
rise to a shift up to 1.5 GeV in the prediction for nth- 

The complete one-loop calculation together with the leading two-loop corrections and 
the other corrections beyond 0(aa s ) have been implemented into the FORTRAN code Feyn- 
Higgs p4 . This code can be linked to existing programs as a subroutine, thus providing an 



accurate calculation of m\ t and m# which can be used for further phenomenological analyses. 

FeynHiggs is available via its WWW page 

|http : / /www-itp . physik . uni-karlsruhe . de/f eynhiggs . 



4 Numerical results for rrih and m# 

4.1 Dependence of the results on the MSSM parameters 

In this subsection we give a detailed discussion of the dependence of m h on the parame- 
ters of the MSSM. For tan f3 we restrict ourselves to two typical values which are favored 



by SUSY-GUT scenarios [27|: tan/5 = 1.6 for the 577(5) scenario and tan/? = 40 for 
the 50(10) scenario. Other parameters are M z = 91.187 GeV, M w = 80.375 GeV, Gf = 
1.16639 10" 5 GeV" 2 , a s {m t ) = 0.1095, m t = 175 GeV, and m b = 4.5 GeV (if not indicated 
differently). The parameter M appearing in the plots is the SU(2) gaugino mass parameter. 
The other gaugino mass parameter, M 1; is fixed via the GUT relation 

M 1 = 5 -^M . (84) 

"J c w 

The scalar top masses and the mixing angle are related to the parameters M^ L , Mf and 
M t LR of the t mass matrix, which reads 

KA 2 _( Ml+ml + cos 2/5 (| - !4,)Mi m t M t LR \ 

JVl i-{ mt M t LR M? r + m\ + § cos 2(3 s 2 w M 2 z J ' [ ^ } 



with 



M t LR = A t - fx cot . (86) 



2 The functional dependence of m t (mt) is known up to 0(a 2 s ) Since m t (mt) enters only at the two- 
loop level, we have incorporated only the one-loop correction to mt(mj), thus neglecting only contributions 
of 0{aa 3 s ) in m,h- 
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In the figures below we have chosen rriq = M^ L = Mf (if not indicated differently). 

Fig. ^ shows the result for obtained from the diagrammatic calculation of the full one- 
loop and leading two-loop contributions. The two contributions beyond 0(otot s ) discussed 
above are shown in separate curves. For comparison the pure one-loop result is also given. 
The results are plotted as a function of M^ R /rriq, where rrtq is fixed to 500 GeV. The two- 
loop contributions give rise to a large reduction of the one-loop result of 10-20 GeV. The two 
corrections beyond 0(aa s ) both increase by up to 2 GeV. A minimum occurs around 
M^ R = GeV which we refer to as 'no mixing'. A maximum in the two-loop result for rrih 
is reached for about M^ R /m, ~ ±2 in the tan (3 =1.6 scenario as well as in the tan (3 = 40 
scenario. This case we refer to as 'maximal mixing'. The position of the maximum is shifted 
compared to its one-loop value of about M^ R /rriq »s ±2.4. The Yukawa correction and 
the insertion of the running top mass have only a negligible effect on the location of the 
maximum. 
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Figure 6: One- and two-loop results for rrih as a function of M^ R /rriq for two values of tan (3 . 
The corrections beyond 0(aa s ) discussed in the text are shown separately. 

Fig. |7] depicts the result for the heavy Higgs-boson mass, tuh, obtained in the same way 
as rrih above. The only difference is that no Yukawa term has been included. In the plot 
we have chosen the small value Ma = 75 GeV, close to the lower experimental bound, since 
only for a light A boson the higher-order corrections give a sizable contribution (see also 
Fig. |J). Here the values for m# obtained for small tan/? are larger than for tan/? = 40. 
The values of M^ R /rriq for which run is maximal depend in this case on tan/3 and the sign 
ofM t LR . 
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Figure 7: One- and two-loop results for ran as a function of M^ R /rriq for two values of tan (3 . 
The running top mass correction discussed in the text is shown separately. 

Both Higgs-boson masses are shown in Fig. |8| for low and high tan (3 and the no-mixing 
and the maximal- mixing case, where the latter case corresponds to the definition according 
to Fig. |6| for the light Higgs boson. For rriu sizable corrections at the one- and two-loop level 
are obtained only for M A < 200 GeV for tan (3 = 1.6 and for M A < 120 GeV for tan (3 = 40. 

More relevant for todays' colliders is the mass of the lighter Higgs boson, rrih, on which 
we will focus in the following discussion. In Fig. |9| is shown in the two scenarios with 
tan/3 = 1.6 and tan/3 = 40 as a function of rriq for no mixing and maximal mixing and 
for Ma = 200, 1000 GeV. The tree-level, the one-loop and the two-loop results with the 
two corrections beyond 0(aa s ) are shown (the values of rriq are such that the corresponding 
t-masses lie within the experimentally allowed region). In all scenarios of Fig. [5] the two- loop 
corrections give rise to a large reduction of the one-loop value of m^. The effect is generally 
larger in the tan/3 =1.6 scenario, and for maximal mixing and large Ma- The inclusion 
of the Yukawa correction and the running top mass leads to a slight shift in rrih towards 
higher values. This effect amounts up to 20% of the two-loop correction. In the tan/3 = 1.6 
scenario with rriq = 1 TeV, rrih reaches about 75 (81) GeV for Ma = 200 (1000) GeV in 
the no-mixing case, and 94 (101) GeV in the maximal-mixing case. For tan/3 = 40 the 
respective values of rrih are 112 GeV in the no-mixing case, and 126 GeV in the maximal- 
mixing case for both values of Ma- The peaks in the plots for Ma = 1 TeV are due to the 
threshold Ma = ro^ + rn~ b2 in the one-loop contribution, originating from the sbottom-loop 
diagram in the A self-energy. 
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Figure 8: m^ and rxiu are shown as a function of Ma for tan/3 = 1.6 and tan (3 = 40. 
The two-loop value for m^ contains both corrections beyond 0(aa s ), whereas for m# only 
the running top mass effect is included. The maximal-mixing scenario corresponds to the 
definition according to the discussion of Fig. ||. 

The dependence of rrih on Ma is depicted in Fig. |1(] in the two scenarios with tan f3 =1.6 
and tan f3 = 40 for no mixing and maximal mixing for uiq = 500, 1000 GeV. In all scenarios 
of Fig. [10] the two-loop corrections give rise to a large reduction of the one-loop value of 
m h . A saturation effect can be observed for Ma ^ 300 (150) GeV in the tan/3 = 1.6 (40) 
scenario. The peaks in the plots for Ma = 350 GeV are due to the threshold Ma = 2 m t in 
the one-loop contribution, originating from the top-loop diagram in the A self-energy. 

Allowing for a splitting between the parameters M^ L , M~ tR in the t mass matrix yields 
maximal values of m^ which are approximately the same as for the case rriq = M^ L = Mr, 
provided that one sets 



Mr =Mr 
1 l L f R 



max{M^,M^} , (87) 



see Fig. 11. However, the location of the maximal Higgs-boson mass, depending on M^ R , is 
shifted towards smaller values, typically by about 40%. The numerical difference in in 
the two splitting scenarios M~ tL /M~ tR = 300/1000 and M~ tL /M~ tR = 1000/300 is small. They 
differ by up to 2 GeV only in the large tan (5 scenario when M^ R > 1000 GeV. 

The variation of with m t is rather strong. The scenarios for no mixing and maximal 
mixing and for tan f3 =1.6 and tan j3 = 40 are shown in Fig. [E| where mt is varied around 
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Figure 9: The mass of the lightest Higgs boson for tan/5 = 1.6 and tan/5 = 40. The 
tree-level, the one-loop and the two-loop results for m h are shown as a function of rriq for 
the no-mixing and the maximal-mixing case, and for Ma = 200, 1000 GeV. 
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Figure 10: The mass of the lightest Higgs boson as a function of M4 for tan/3 = 1.6 and 
tan/5 = 40. The tree-level, the one-loop and the two-loop results for m h are shown in the 
no-mixing and the maximal-mixing case, and for rriq = 500, 1000 GeV. 
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Figure 11: The mass of the lightest Higgs boson as a function of M t for differently split 
values of the soft SUSY breaking terms. The curves in the plots correspond to the values 
M- tL /M- tR = 1000/1000, 300/1000 and 1000/300. 
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the central value of m t = 175 GeV by ±10 GeV. The variation of rrih is stronger for low 
tan f3 and larger trig-, in the tan {3 =1.6 scenario rrih varies by more than 10 GeV and about 
20 GeV for no-mixing and maximal-mixing, respectively. In the tan/5 = 40 scenario the 
respective values are less than 10 GeV and about 15 GeV. 

Varying tan/? around the value tan/5 = 1.6 has a relatively large impact on rrih (higher 
values for rrih are obtained for larger tan (3 ), while the effect of varying tan f3 around tan (3 = 
40 is marginal. This is shown in Fig. |TJ for M A = 200, 1000 GeV, m- q = 200, 1000 GeV for 
the no-mixing and the maximal-mixing scenario. For tan (3 > 15 the variation is less than 
1 GeVQ. 



In Fig. O rrih is shown as a function of M, the soft SUSY breaking parameter in the 



chargino and neutralino sector (see Sec. 84). In our calculation M enters only in the one- 
loop self-energies. The variation is less than 4 GeV for the whole M parameter space. For 
increasing M the result for rrih decreases in general. 

The dependence of rrih on the Higgs mixing parameter, is depicted in Fig. |T5| The 
parameter /i enters via the non-diagonal Higgs-squark coupling at one- and two-loop order 
and via the chargino and neutralino sector in the one-loop self-energies. It should be noted 



that for the plots in Fig. [L5| we have set m& = GeV, thus suppressing the contribution of the 
b — 6-sector. The reason is that for large \i and for large tan (3 some Higgs-sbottom couplings 
can become rather large, which makes the perturbative calculation questionable in this case. 



The variation of rrih with \i in Fig. [15| is relatively weak, not exceeding 3 GeV. A maximum 



(for the choice M = 400 GeV) for m h lies between /i = -200 GeV and fi = -100 GeV. For 



decreasing M the maximum is reached for slightly smaller values of /z, see also Sec. |4~2 . 



Finally we show the dependence on the gluino mass, rrig, which enters only at the two- 
loop level. Fig. [16] depicts rrih as a function of m g in the scenarios with rrig = 500, 1000 GeV 
for tan (3 = 1.6 and tan/? = 40 in the no-mixing and the maximal-mixing case. Small 
variations below 1 GeV occur in the no-mixing scenario, while change in rrih up to 4 GeV 
arises in the maximal-mixing scenario, rrih reaches a maximum at about m g »s 0.8 rrig. Since 
the parameter m g is absent in the RG approach, a variation of rrih with m g can directly be 
seen as a deviation of the diagrammatic result from the RG result, see Sec. |1.3|. 



As pointed out in Ref. it is desirable to express the predictions for the observable 
rrih in terms of other physical observables. This provides the possibility to directly compare 
results obtained by different approaches making use of different renormalization schemes. 
Therefore we show in Fig. [H] the dependence of rrih on the parameters , mi 2 and 6+, 
which, since we are working in the on-shell scheme, directly correspond to the physical 
ones. We show rrih as a function of m t ~ 2 for Am t - = GeV and 6~ t = (no mixing) and 
for Am; = 340 GeV and Q~ t = —it/A (maximal mixing), where Am^ = m t ~ 2 — m~ tl . The 
choice of Ami ~ 340 GeV corresponds to M^ R /rriq « 2 in terms of the soft SUSY breaking 

3 A non-negligible effect can arise for large tan /3 if (x is also large. This is briefly discussed below in the 
context of the /^-dependence of mj,. 



23 



140 
120 
100 
® 80 

a 

E 60 
40 
20 



m~ = 500 GeV, m~ = 500 GeV, M A = 500 GeV, tanp = 1.6 



— tree level 

1 -loop, no mixing 

2-loop, no mixing, incl. all ret. 

1 -loop, max. mixing 

2-loop, max. mixing, incl. all ref. 





165 167 169 171 173 175 177 179 181 183 185 
mJGeV] 



160 : 
140 - 
120 

1 100 ; 

E 80 -m~ = 500 GeV,m~ = 500GeV,M A = 500GeV,tanp = 40 : 

60 - — tree level 

1 -loop, no mixing 

2-loop, no mixing, incl. all ret. 
40 - 1-loop, max. mixing 

2-loop, max. mixing, incl. all ref. 

2Q r 

165 167 169 171 173 175 177 179 181 183 185 
mJGeV] 




— tree level 

1-loop, no mixing 

2-loop, no mixing, incl. all ret. 

1-loop, max. mixing 

2-loop, max. mixing, incl. all ref. 



165 167 169 171 173 175 177 179 181 183 185 
mJGeV] 



160 : 
140 - 
120 



\ 100 



E 80 - m~ = 1 000 GeV, m- = 500 GeV, M A = 500 GeV, tanp = 40 



40 



20 



— tree level 

1-loop, no mixing 

2-loop, no mixing, incl. all ref. 

1-loop, max. mixing 

2-loop, max. mixing, incl. all ref. 



165 167 169 171 173 175 177 179 181 183 185 
mJGeV] 



Figure 12: The mass of the lightest Higgs boson is shown as a function of m t : the 
scenarios of no mixing and maximal mixing are depicted for low and high tan (5 , and 
fx = -200 GeV, M = 400 GeV. 
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Figure 13: The mass of the lightest Higgs boson is shown as a function of tan/3 . m h is 
plotted for the scenarios with M A = 200, 500 GeV and mq = 500, 1000 GeV in the no-mixing 
and the maximal-mixing case. 
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Figure 14: The mass of the lightest Higgs boson is shown as a function of M, the soft SUSY 
breaking parameter of the chargino and neutralino sector. 



parameters. The two-loop results shown here contain also the corrections beyond 0(aa s ). 
In these plots we have furthermore imposed the p -parameter constraint: We have required 
that the contribution of the third generation of scalar quarks to the p-parameter ^] 
does not exceed the value of 1.3 • 10~ 3 , which corresponds approximately to the resolution of 
Ap when it is determined from experimental data [ 2"8fl . For tan/3 = 1.6 reaches about 
76 (82) GeV for M A = 200 (1000) GeV,m f ~ 2 = 1 TeV and no mixing in the t-sector. In 
the maximal-mixing case the reached values of rrih are 94 (101) GeV. In the tan/3 = 40 
scenario, rrih reaches 114 (127) GeV in the no-mixing (maximal-mixing) case for both values 
of Ma- The peaks in the plots for M A = 1 TeV and maximal mixing in the t-sector around 
mi 2 = 660 GeV are due to the threshold Ma = + m^ 2 in the one-loop contribution, 
originating from the stop-loop diagram in the A self-energy. 

4.2 Upper bound for rrih as a function of tan/3 

Since, as shown in Fig. [13|, smaller values for rrih are obtained for small tan f3 , this part of 
the parameter space can to a large extend be covered at todays' colliders. The discovery 
limit for rrih at LEP2 is expected to be slightly above 100 GeV |[29|| . In this context it is of 



special interest to know the maximally possible value for rrih as a function of tan (3 in the 
MSSM. To this end we have performed a parameter scan, varying rrig, M, p, M A and M t LR for 
three values of m t and fixed values of rrig and tan (3 . The maximal values for rrih, including 
also the Yukawa correction and the contribution from the running top mass, were reached 
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Figure 15: The mass of the lightest Higgs boson is shown as a function of /x, the Higgs 
mixing parameter. The mass of the bottom quark, mj, is set to zero in order to avoid 
unnaturally large one-loop effects. 

(in the case m t = 175 GeV) forf] 

m- g pa 0.8 rriq 
M » GeV 
fi « GeV 



M A 



M t LR Pi 2m g (88) 
in all scenarios. 

The value for Ma in fl88p needs some further explanation: as one can see in Fig. [TT], due 
to a one-loop threshold effect the value of m^ can become very large for Ma = 2m t . For 
large tan f3 this threshold effect results in a bigger rrih value for Ma in the region around 2 mt 
than for larger values of Ma (with Ma < 1500 GeV, where we stopped our scan.) Of course 
the exact value Ma = 2 m t = 2 ■ 175 GeV = 350 GeV would be a very specific choice, giving 
a wrong impression of the possible size of m^. (Exactly at the threshold also finite width 
effects for the A boson would have to be taken into account.) Therefore we have chosen 
the value Ma = 360 GeV which is not directly at the threshold, thus giving a more realistic 



800 GeV for small tan/? 
360 GeV for large tan (5 



4 Due to threshold effects very high values for rrih can occur. Since this is regarded as an accidental effect, 
these isolated points of parameter space are not considered here. 
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Figure 16: The mass of the lightest Higgs boson as a function of m~ g for a common value 
of M A = 500 GeV, rriq = 500, 1000 GeV for the no-mixing and the maximal-mixing case and 
for low and high tan (5 . 
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Figure 17: The mass of the lightest Higgs boson in terms of the physical parameters 
m h > m h an d ®b where Am ( - = mi 2 — . The scenarios Am^ = GeV, (% = (no mixing) 
and Am; = 340 GeV, 6i = — tt/4 (maximal mixing) are shown. 
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impression about the maximally possible values for m^. The choice M = \i = GeV is 
experimentally excluded. We nevertheless use these values since the difference in rrih to the 
case with experimentally not excluded M and /i is very small, typically below 0.5 GeV. 

In Fig. we show the maximal Higgs-boson mass value, including also the corrections be- 
yond 0(aa s ), as a function of tan/3 ; the other parameters are chosen according to eqs. fl8*B|) . 
For the top-quark mass the most recent experimental value m t = 173.8 GeV |Kj is chosen 
and, since rrih grows with increasing m t , the experimental value plus one and plus two stan- 
dard deviations (m t = 178.8, 183.8 GeV)Q. The common squark mass parameter is chosen 
to be rriq = 1000 GeV as a high, and rriq = 2000 GeV as a very high value. On the left side 
of Fig. [18] we show the full tan/3 range (tan/3 < 50), whereas on the right side we focus on 
the range especially interesting for LEP2 and the upgraded Tevatron (tan/3 < 5). 

In the tan/3 < 5 plot we have chosen Ma = 800 GeV. In the tan/3 < 50 plot, however, 
we have chosen Ma = 800 GeV only for tan/3 < 4; for larger values we have switched to 
Ma = 360 GeV. For the value tan (3 = 4 one gets about the same maximal value for rrih for 
both choices of Ma- 
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Table 1: Maximal values for rrih for different choices of m t ,mq and tan/3 . All other para- 



meters have been chosen according to eqs. (88). (All masses are in GeV. 



In the interesting region around tan/3 = 1.6 the covered region of the tan (3 -parameter 
space depends strongly on the maximally accessible energy of todays' colliders, see Tab. [l|. 
For an exclusion limit of rrih > 107 GeV, for instance, LEP2 covers tan/3 < 1.6 completely 
only if m t is constrained to its present 1 a limit. On the other hand, taking a very conservative 
point of view and choosing m t at the 2a bound, no limit on tan/3 can be set, even for 
rriq = 1000 GeV. 



One should keep in mind, however, that the Higgs-boson masses depicted in Fig 
are the maximally possible upper values, i.e. for smaller mixing in the t-sector the region 
tan/3 < 1.6 can be covered by LEP2 for all other sets of parameters. One can also see that 
a precise measurement of m t is decisive in order to set stringent bounds on tan/3 in the 
MSSM. 



One should note, however, that the highest value for mt is disfavored in the MSSM by internal consis- 
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In conclusion, our results confirm that for the scenario with tan/3 = 1.6 the parameter 
space of the MSSM can be covered to a very large extent. Only for maximal mixing, very 
large soft SUSY breaking parameters in the i-sector and m t at its upper (1 — 2)a limit the 
light Higgs boson can escape the detection at LEP2. For increasing tan f3 , however, the 
parameter space in which the Higgs boson is not accessible at LEP2 increases rapidly. 

Concerning the large tan f3 region, LEP2 and the upgraded Tevatron can probe only the 
region of no mixing in the t-sector. The LHC and a future linear e + e~-collider are needed 
in order to test the parameter space with large t-mixing. 

In an analogous way we have also analyzed the maximal value for as a function of the 



physical parameters, see Fig. [19] and Tab. ||. We have chosen m ( - 2 = 1000, 2000 GeV, Am ( - 



340 GeV and 9$ = — 7r/4. The other MSSM parameters are chosen according to eqs. 
Fig. shows the result for the maximal range of tan/5 and for the interesting range for 
LEP2, tan/3 < 5. The results in Fig. pi] and Tab. are slightly lower than the values 
obtained with the unphysical input parameters. This is due to the fact that the values 
obtained for the squark masses in the first scenario (for all other parameters chosen to be 
equal) are always larger than for the latter case with physical input parameters. The analysis 
of the upper bound of rrih, however, can be taken over directly from the case with unphysical 
input parameters. 
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Table 2: Maximal values for for different choices of m t ,mf and tan/3 . All other para- 
meters have been chosen according to eqs. (|88|). (All masses are in GeV.) 



4.3 Numerical comparison with the RG approach 

We now turn to the comparison of our diagrammatic results with the predictions obtained 
via RG methods. For this comparison we made use of the FORTRAN code corresponding 
to Ref. |§, except for the one-loop results in Figs. |2(| and [H], where we used the code 
described in Ref. |nj[J 

6 The RG results of Ref. § and Ref. go| agree within about 2 GeV with each other. 
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Figure 18: The maximally possible value for uih, including all refinement terms, as a 
function of tan/? , depending on the unphysical parameters nig and M^ R . The other MSSM 
parameters have been chosen according to eqs. flSSp . 
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Figure 19: The maximally possible value for m^, including all refinement terms, as a 
function of tan (3 , depending on the physical parameters m 4 ~ 2 , Am t - and %. The other MSSM 
parameters have been chosen according to eqs. (jSSp . 
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We begin with the case of large values of M4, for which the RG approach is most easily 
applicable and is expected to work most accurately. In order to study different contributions 
separately, we have first compared the diagrammatic one-loop on-shell result with the 
one-loop leading log result (without renormalization group improvement) given in Ref. |10| . 



Since the available code uses the running top mass m t = m t (m t ) m 167.3 GeV we have also 
used this top mass for the full diagrammatic one-loop calculation. In Fig. 2C the lightest 
Higgs-boson mass is shown in the no-mixing scenario, i.e. M^ R = GeV, whereas in Fig. ^T] 
rrih is shown for increasing mixing in the t-sector. We found very good agreement, typically 
within 1 GeV for both mixing cases and low and high tan (3 . Only for very small values of 
rriq a deviation up to 2 GeV arises. For values of M4 below 100 GeV (which are not shown 
here) and large mixing in the t-sector deviations of about 5 GeV occur. 
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Figure 20: Comparison between the one-loop Feynman-diagrammatic calculations and the 
results obtained by renormalization group methods ]TU|] . The mass of the lightest Higgs 
boson is shown for the two scenarios with tan/3 = 1.6 and tan/3 = 40 as a function of 
M t LR /mq for M A = 200, 1000 GeV. 

In the next step of comparison we analyzed the no-mixing case at the two-loop level: 
we have compared our diagrammatic result for the no-mixing case, including the Yukawa 
correction and the running top mass effect, with the RG results obtained in Ref. f|. We 
have adopted the scale M (the SU{2) gaugino mass parameter) as M = m^, in order to 



"</• 



treat it in the same way as it has been done in the RG approach. As can be seen in Fig. 
after the inclusion of the corrections beyond 0(aa s ) the diagrammatic result for the no- 
mixing case agrees very well with the RG result. For the scenario with M4 = 1000 GeV the 
deviation between the results exceeds 2 GeV only for tan j3 



1.6 and nig < 150 GeV. For 



34 





160 




140 




120 




100 












o 


80 


e" 






60 




40 




20 








-3.0 



m- = 1000 GeV, m g ~ = 500 GeV, n = -m- 
M A = 200 GeV, m, = 1 66.5 GeV 




RG 1-loop, tanp= 1.6 
full 1 -loop, tanfj = 1 .6 
RG 1-loop, tanp = 40 
full 1-loop, tan|5 = 40 



-2.0 



-1.0 



0.0 

Mf/m- 



1.0 



2.0 



3.0 





160 




140 




120 




100 


> 









o 


80 


E 






60 




40 




20 








-3.0 



m- = 1000 GeV, m- = 500 GeV, u. = -m- 

q g > r q 



M A = 1000 GeV, m 



166.5 GeV 




RG 1-loop, tanp= 1.6 
full 1-loop, tanp = 1.6 
RG 1-loop, tanp = 40 
full 1-loop, tan|3 = 40 



-2.0 



-1.0 



0.0 

r</m~ 



1.0 



2.0 



3.0 



Figure 21: Comparison between the one-loop Feynman-diagrammatic calculations and the 
results obtained by renormalization group methods flOf . The mass of the lightest Higgs 
boson is shown for the two scenarios with tan/? = 1.6 and tan/3 = 40 as a function of 
M t LR /mg for M A = 200, 1000 GeV. 
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Ma = 200 GeV the deviation is in general slightly larger than for Ma = 1000 GeV, but does 
not exceed 2 GeV. 

The RG results do not contain the gluino mass as a parameter. Hence, varying m§, which 



has been discussed in Sec. [4J], gives rise to an extra deviation. In the no-mixing case this 
extra deviation does not exceed 1 GeV. Varying the other parameters /i and M in general 
does not lead to a sizable effect in the comparison with the corresponding RG results (as 
long as M^ R is taken as input and a variation of /j, does not affect the t-mixing.) 

Finally we consider the situation where mixing in the t sector is taken into account. 
In Fig. |23| our diagrammatic result, including the Yukawa correction and the running top 
mass effect, is compared with the RG results as a function of M^ R /rriq for the cases 
tan/3 = 1.6 and tan/? = 40, and for m q = 500,1000 GeV and M A = 200 GeV. The 



Ma = 1000 GeV scenario is depicted in Fig. ^ for the same set of parameters. The point 
M^ R /nig = corresponds to the plots shown in Fig. |22], except that the parameter // is 
set to fi — —trig here. For larger t-mixing, sizable deviations between the diagrammatic 
and the RG results occur. They can reach 5 GeV for moderate mixing and become very 
large for \M^ R /rriq\ > 2.5. As already mentioned above, the maximal value for rrih in 
the diagrammatic approach is reached for M^ R /mq ~ ±2, whereas the RG results have a 
maximum at M t LR /m q « ±2.4, i.e. at the one-loop value. This holds for all combinations 
of tan/3 , rriq and Ma- In the case of positive M^ R , the maximal values for reached in 
the diagrammatic calculation are up to 5 (3) GeV larger than the ones of the RG method 
for tan/5 =1.6 (40). The dependence on M f LR is asymmetric; for negative M t about the 
same maximal values are reached in the two approaches. 

The diagrammatic result varies with rrig as shown in Fig. 16. In the case of mixing in 



the t-sector this leads in general to a larger effect than in the no-mixing case and shifts the 
diagrammatic result relative to the RG result within ±2 GeV. 

Up to now we have compared the results of our diagrammatic on-shell calculation and 
the RG methods in terms of the (unphysical) soft SUSY breaking parameters of the t mass 
matrix M$ L , M^ R and M^ R , since the available numerical codes for the RG results || 1C] 



are given in terms of these parameters. However, since the two approaches rely on different 
renormalization schemes, the meaning of these non-observable parameters is not precisely 
the same in the two approaches starting from two-loop order. Indeed we have checked 
that assuming fixed values for the physical parameters m^, m 4 ~ 2 , and 9i and deriving the 
corresponding values of the parameters Mj , M- tR and M^ R in the on-shell scheme as well 
as in the MS scheme, sizable differences occur between the values of the mixing parameter 
M^ R in the two schemes. On the other hand the parameters Mf L , M^ R are approximately 
equal in both schemes. Thus, part of the different shape of the curves in Fig. |23| and Fig. p4| 
may be attributed to a different meaning of the parameter M^ R in the on-shell scheme and 
in the RG calculation. 

In order to avoid this problem in comparing results obtained by different approaches 
making use of different renormalization schemes, we find it preferable to compare predictions 
for physical observables in terms of other observables (instead of unphysical parameters). 
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Figure 22: Comparison between the Feynman-diagrammatic calculations and the results 
obtained by renormalization group methods |J. The mass of the lightest Higgs boson is 
shown for the scenarios with tan/3 = 1.6 and tan/3 = 40, Ma = 200, 1000 GeV for the case 
of vanishing mixing in the t-sector. 
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Figure 23: Comparison between the Feynman-diagrammatic calculations and the results 
obtained by renormalization group methods |J. The mass of the lightest Higgs boson is 
shown for the two scenarios with tan/3 = 1.6 and tan/5 = 40 as a function of M^ R /m,q for 
nig = 200, 1000 GeV and M A = 200 GeV. 
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Figure 24: Comparison between the Feynman-diagrammatic calculations and the results 
obtained by renormalization group methods |J. The mass of the lightest Higgs boson is 
shown for the two scenarios with tan/3 = 1.6 and tan/5 = 40 as a function of M^ R /m,q for 
nig = 200, 1000 GeV and M A = 1000 GeV. 
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Therefore we switch from the set of unphysical parameters to a set of physical parameters: 



M h ,M iR ,M t LR - m 6 ,,Am t -(=m fe -m fl ),0 £ . (89) 



In Fig. |25]we compare the results for the lightest Higgs-boson mass, obtained by the Feynman- 
diagrammatic method and by the RG method, in terms of this new set of parameters: 
is shown as a function of m: with the mass difference Am t - = mr — m £ and the mixing 
angle df as further input parameters. In the context of the RG approach the running t- 
masses, derived from the t mass matrix, are considered as an approximation for the physical 
masses. In our approach, on the other hand, since we are working in the on-shell scheme, 
the t-masses and the mixing angle directly correspond to physical parameters. In Fig. |25] we 
have furthermore implemented the same Ap constraints on the range of the third generation 
scalar quark masses as in Fig. |17]. 



Similarly to the comparison shown in Fig. and |23, very good agreement is found 



in Fig. |25| between the results of the two approaches in the case of vanishing t-mixing. 
The deviation is typically less than 1 GeV and never exceeds 2 GeV. Using the physical 
parameters as input, the maximal-mixing scenario is realized by setting 9i = —tt/A and 
Affl ( - « 340 GeV (i.e. the t-masses obtained for M^ R /rriq « 2 have a mass difference of 
about 340 GeV.) In this scenario again (as in Figs. ^3]and|24|) the diagrammatic result yields 
values for rrih which are higher by about 5 GeV. The peaks in the plots for Ma = 1 TeV 
and maximal mixing in the i-sector around = 660 GeV are again due to the threshold 
Ma = Wj + m ( - 2 in the one-loop contribution, originating from the stop-loop diagram in the 
A self-energy. 



5 Conclusions 

Using the Feynman diagrammatic method we have calculated the leading 0(aa s ) corrections 
to the masses of the neutral CP-even Higgs bosons in the MSSM. The two-loop result has 
been implemented into the prediction based on the complete diagrammatic one-loop on-shell 
result. Two further corrections beyond 0(aa s ) have been added in order to incorporate 
leading electroweak two-loop and higher-order QCD contributions. The results have been 
obtained using the on-shell scheme, which means a renormalization of all sectors of the MSSM 
at one-loop order and of the Higgs-boson sector at two-loop order. In our two-loop calculation 
we have imposed no restrictions on the parameters of the Higgs and scalar top sector of the 
model. Thus the results are valid for arbitrary values of the relevant MSSM parameters. 
The complete result has been implemented into the FORTRAN program FeynHiggs |24| which 
is available via its WWW page 

trttp : / /www-itp . physik . uni-karlsruhe . de/f eynhiggsj . 

In this way we provide the at present most precise prediction for and m# based on 
Feynman-diagrammatic calculations. 

The two- loop corrections lead to a large reduction of the one- loop on-shell result. We have 
performed a detailed analysis of the dependence of on the various MSSM parameters. 
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Figure 25: Comparison between the Feynman-diagrammatic calculations and the results 
obtained by renormalization group methods || in terms of physical parameters. The mass 
of the lightest Higgs boson is shown for the two scenarios with tan/3 =1.6 and tan/3 = 40 
and for Ma = 200, 1000 GeV as a function of the heavier physical t mass mi . For the curves 
with 0j = a mass difference Am ( - = GeV is taken, whereas for 9f = —ir/4 we choose 
Am f - = 340 GeV, for which the maximal Higgs-boson masses are achieved. 
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Concerning the scalar top sector the analysis has been carried out in terms of the (unphysical) 
soft SUSY breaking parameters M^ L , Mi R and M^ R as well as in terms of the physical 
parameters mi 2 , Am ( - = mf 2 — m~ tx and B~ t . 

A scan over the parameters fj,, M,rrig, Ma and M^ R has been performed in order to 
determine the maximally possible value for rrih as a function of tan j3 . Our results show 
that for the scenario with tan/3 = 1.6 the parameter space of the MSSM can be covered 
almost completely. Only for maximal mixing, very large soft SUSY breaking parameters 
in the t-sector and m t at its upper experimental limit the light Higgs boson can escape 
the detection at LEP2 in this scenario. Concerning the large tan j3 region, LEP2 and the 
upgraded Tevatron can probe only the region of no mixing in the t-sector. 

We have compared our results, obtained by a Feynman diagrammatic calculation (where 
also the corrections beyond 0(aa s ) have been included), with the results obtained via RG 
methods. Concerning the one-loop contributions we find very good agreement between 
these two approaches. The same is valid for the two-loop corrections in the case of vanishing 
mixing in the t-sector. On the other hand, in the case of non-vanishing mixing sizable 
deviations between the two approaches occur. For moderate mixing they reach up to 5 GeV, 
for \M t LR /mq\ > 2.5 they can be very large. In the diagrammatic approach the maximal 
value for m h is reached for M^ R /rriq ps ±2, whereas the RG results have a maximum at 
M^ R /rriq ps ±2.4, i.e. at the one-loop value. This holds for all combinations of tan/5 ,m 9 - 
and Ma- The fact that the parameter m g is absent in the RG results can give rise to an 
additional deviation between the two approaches of about ±2 GeV. 

We have furthermore discussed the issue of how results obtained via different approaches 
using different renormalization schemes can be readily compared to each other also when 
corrections beyond one-loop order are incorporated. For this purpose it is adequate to 
express the prediction for the Higgs-boson masses in terms of other physical observables, i.e. 
the physical masses and mixing angles of the model. 

Accordingly, we have compared the results obtained by our diagrammatic two-loop cal- 
culation with those obtained by RG methods in terms of the physical observables m t ~ 2 , 
Ami = mi — mi and % As for the comparison in terms of the unphysical parameters, 
we have found good agreement for the case of vanishing mixing in the t-sector. For large 
splitting between the t-masses, however, the Higgs-boson masses obtained by the Feynman 
diagrammatic calculation are about 5 GeV larger than the ones calculated in the RG ap- 
proach. 
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